Delayed referral for liver transplant evaluation worsens outcomes in chronic liver disease patients requiring inpatient transplant evaluation.

BACKGROUND: Indications to refer patients with cirrhosis for liver transplant evaluation (LTE) include hepatic decompensation or a model for end stage liver disease (MELD-Na) score ≥ 15. Few studies have evaluated how delaying referral beyond these criteria affects patient outcomes. AIM: To evaluate clinical characteristics of patients undergoing inpatient LTE and to assess the effects of delayed LTE on patient outcomes (death, transplantation). METHODS: This is a single center retrospective cohort study assessing all patients undergoing inpatient LTE (n = 159) at a large quaternary care and liver transplant center between 10/23/2017-7/31/2021. Delayed referral was defined as having prior indication (decompensation, MELD-Na ≥ 15) for LTE without referral. Early referral was defined as referrals made within 3 mo of having an indication based on practice guidelines. Logistic regression and Cox Hazard Regression were used to evaluate the relationship between delayed referral and patient outcomes. RESULTS: Many patients who require expedited inpatient LTE had delayed referrals. Misconceptions regarding transplant candidacy were a leading cause of delayed referral. Ultimately, delayed referrals negatively affected overall patient outcome and an independent predictor of both death and not receiving a transplant. Delayed referral was associated with a 2.5 hazard risk of death. CONCLUSION: Beyond initial access to an liver transplant (LT) center, delaying LTE increases risk of death and reduces risk of LT in patients with chronic liver disease. There is substantial opportunity to increase the percentage of patients undergoing LTE when first clinically indicated. It is crucial for providers to remain informed about the latest guidelines on liver transplant candidacy and the transplant referral process.

Characterization of triatomine bloodmeal sources using direct Sanger sequencing and amplicon deep sequencing methods.

Knowledge of host associations of blood-feeding vectors may afford insights into managing disease systems and protecting public health. However, the ability of methods to distinguish bloodmeal sources varies widely. We used two methods-Sanger sequencing and amplicon deep sequencing-to target a 228 bp region of the vertebrate Cytochrome b gene and determine hosts fed upon by triatomines (n = 115) collected primarily in Texas, USA. Direct Sanger sequencing of PCR amplicons was successful for 36 samples (31%). Sanger sequencing revealed 15 distinct host species, which included humans, domestic animals (Canis lupus familiaris, Ovis aries, Gallus gallus, Bos taurus, Felis catus, and Capra hircus), wildlife (Rattus rattus, Incilius nebulifer, Sciurus carolinensis, Sciurus niger, and Odocoileus virginianus), and captive animals (Panthera tigris, Colobus spp., and Chelonoidis carbonaria). Samples sequenced by the Sanger method were also subjected to Illumina MiSeq amplicon deep sequencing. The amplicon deep sequencing results (average of 302,080 usable reads per sample) replicated the host community revealed using Sanger sequencing, and detected additional hosts in five triatomines (13.9%), including two additional blood sources (Procyon lotor and Bassariscus astutus). Up to four bloodmeal sources were detected in a single triatomine (I. nebulifer, Homo sapiens, C. lupus familiaris, and S. carolinensis). Enhanced understanding of vector-host-parasite networks may allow for integrated vector management programs focusing on highly-utilized and highly-infected host species.

Amplicon Sequencing as a Potential Surveillance Tool for Complexity of Infection and Drug Resistance Markers in Plasmodium falciparum Asymptomatic Infections.

BACKGROUND: Genotyping Plasmodium falciparum subpopulations in malaria infections is an important aspect of malaria molecular epidemiology to understand within-host diversity and the frequency of drug resistance markers. METHODS: We characterized P. falciparum genetic diversity in asymptomatic infections and subsequent first febrile infections using amplicon sequencing (AmpSeq) of ama1 in Coastal Kenya. We also examined temporal changes in haplotype frequencies of mdr1, a drug-resistant marker. RESULTS: We found >60% of the infections were polyclonal (complexity of infection [COI] >1) and there was a reduction in COI over time. Asymptomatic infections had a significantly higher mean COI than febrile infections based on ama1 sequences (2.7 [95% confidence interval {CI}, 2.65-2.77] vs 2.22 [95% CI, 2.17-2.29], respectively). Moreover, an analysis of 30 paired asymptomatic and first febrile infections revealed that many first febrile infections (91%) were due to the presence of new ama1 haplotypes. The mdr1-YY haplotype, associated with chloroquine and amodiaquine resistance, decreased over time, while the NY (wild type) and the NF (modulates response to lumefantrine) haplotypes increased. CONCLUSIONS: This study emphasizes the utility of AmpSeq in characterizing parasite diversity as it can determine relative proportions of clones and detect minority clones. The usefulness of AmpSeq in antimalarial drug resistance surveillance is also highlighted.

Sensitive, Highly Multiplexed Sequencing of Microhaplotypes From the Plasmodium falciparum Heterozygome.

BACKGROUND: Targeted next-generation sequencing offers the potential for consistent, deep coverage of information-rich genomic regions to characterize polyclonal Plasmodium falciparum infections. However, methods to identify and sequence these genomic regions are currently limited. METHODS: A bioinformatic pipeline and multiplex methods were developed to identify and simultaneously sequence 100 targets and applied to dried blood spot (DBS) controls and field isolates from Mozambique. For comparison, whole-genome sequencing data were generated for the same controls. RESULTS: Using publicly available genomes, 4465 high-diversity genomic regions suited for targeted sequencing were identified, representing the P. falciparum heterozygome. For this study, 93 microhaplotypes with high diversity (median expected heterozygosity = 0.7) were selected along with 7 drug resistance loci. The sequencing method achieved very high coverage (median 99%), specificity (99.8%), and sensitivity (90% for haplotypes with 5% within sample frequency in dried blood spots with 100 parasites/µL). In silico analyses revealed that microhaplotypes provided much higher resolution to discriminate related from unrelated polyclonal infections than biallelic single-nucleotide polymorphism barcodes. CONCLUSIONS: The bioinformatic and laboratory methods outlined here provide a flexible tool for efficient, low-cost, high-throughput interrogation of the P. falciparum genome, and can be tailored to simultaneously address multiple questions of interest in various epidemiological settings.

Plasmodium falciparum is evolving to escape malaria rapid diagnostic tests in Ethiopia.

In Africa, most rapid diagnostic tests (RDTs) for falciparum malaria recognize histidine-rich protein 2 antigen. Plasmodium falciparum parasites lacking histidine-rich protein 2 (pfhrp2) and 3 (pfhrp3) genes escape detection by these RDTs, but it is not known whether these deletions confer sufficient selective advantage to drive rapid population expansion. By studying blood samples from a cohort of 12,572 participants enroled in a prospective, cross-sectional survey along Ethiopia's borders with Eritrea, Sudan and South Sudan using RDTs, PCR, an ultrasensitive bead-based immunoassay for antigen detection and next-generation sequencing, we estimate that histidine-rich protein 2-based RDTs would miss 9.7% (95% confidence interval 8.5-11.1) of P. falciparum malaria cases owing to pfhrp2 deletion. We applied a molecular inversion probe-targeted deep sequencing approach to identify distinct subtelomeric deletion patterns and well-established pfhrp3 deletions and to uncover recent expansion of a singular pfhrp2 deletion in all regions sampled. We propose a model in which pfhrp3 deletions have arisen independently multiple times, followed by strong positive selection for pfhrp2 deletion owing to RDT-based test-and-treatment. Existing diagnostic strategies need to be urgently reconsidered in Ethiopia, and improved surveillance for pfhrp2 deletion is needed throughout the Horn of Africa.

Plasmodium falciparum DHFR and DHPS Mutations Are Associated With HIV-1 Co-Infection and a Novel DHPS Mutation I504T Is Identified in Western Kenya.

Antifolate resistance is significant in Kenya and presumed to result from extensive use and cross-resistance between antifolate antimalarials and antibiotics, including cotrimoxazole/Bactrim used for HIV-1 chemotherapy. However, little is known about antifolate-resistant malaria in the context of newly diagnosed HIV-1 co-infection prior to administration of HIV-1 chemotherapy. Blood samples from a cross-sectional study of asymptomatic adult Kenyans enrolled during voluntary HIV testing were analyzed by PCR for Plasmodium spp. More than 95% of volunteers with identifiable parasite species (132 HIV-1 co-infected) were infected with Plasmodium falciparum alone or P. falciparum with Plasmodium ovale and/or Plasmodium malariae. Deep sequencing was used to screen for mutations in P. falciparum dihydrofolate reductase (dhfr) (N51I, C59R, S108N, I164L) and dihydropteroate synthase (dhps) (S436H, A437G, K540E, A581G) from 1133 volunteers. Individual mutations in DHPS but not DHFR correlated with HIV-1 status. DHFR haplotype diversity was significantly different among volunteers by gender and HIV-1 status. DHPS haplotype diversity by HIV-1 status was significantly different between volunteers paired by age and gender, indicating that patterns of resistance were independent of these variables. Molecular simulations for a novel DHPS mutation (I504T) suggested that the mutated protein has increased affinity for the endogenous ligand DHPPP and decreased affinity for drug binding. A sub-group of monoclonal infections revealed that age and parasitemia were not correlated and enabled identification of a rare septuple-mutant haplotype (IRNL-HGEA). In our study, adult Kenyans newly diagnosed with HIV-1 infection were predominantly infected with moderately resistant P. falciparum, with patterns of infecting parasite genotypes significantly associated with HIV-1 status. Together with the discovery of DHPS I504T, these data indicate that antifolate resistance continues to evolve in Kenya. Further, they highlight the need to understand the effects of associated mutations on both fitness and resistance of P. falciparum in the context of HIV-1 co-infection to better inform treatment for asymptomatic malaria.

The impact of antimalarial resistance on the genetic structure of Plasmodium falciparum in the DRC.

The Democratic Republic of the Congo (DRC) harbors 11% of global malaria cases, yet little is known about the spatial and genetic structure of the parasite population in that country. We sequence 2537 Plasmodium falciparum infections, including a nationally representative population sample from DRC and samples from surrounding countries, using molecular inversion probes - a high-throughput genotyping tool. We identify an east-west divide in haplotypes known to confer resistance to chloroquine and sulfadoxine-pyrimethamine. Furthermore, we identify highly related parasites over large geographic distances, indicative of gene flow and migration. Our results are consistent with a background of isolation by distance combined with the effects of selection for antimalarial drug resistance. This study provides a high-resolution view of parasite genetic structure across a large country in Africa and provides a baseline to study how implementation programs may impact parasite populations.

Antimalarial Drug Resistance Profiling of Plasmodium falciparum Infections in Ghana Using Molecular Inversion Probes and Next-Generation Sequencing.

A key drawback to monitoring the emergence and spread of antimalarial drug resistance in sub-Saharan Africa is early detection and containment. Next-generation sequencing methods offer the resolution, sensitivity, and scale required to fill this gap by surveilling for molecular markers of drug resistance. We performed targeted sequencing using molecular inversion probes to interrogate five Plasmodium falciparum genes (pfcrt, pfmdr1, pfdhps, pfdhfr, and pfk13) implicated in chloroquine, sulfadoxine-pyrimethamine (SP), and artemisinin resistance in two sites in Ghana. A total of 803 dried blood spots from children aged between 6 months and 14 years presenting with uncomplicated P. falciparum malaria at the Begoro District Hospital in Begoro and the Ewim Polyclinic in Cape Coast, Ghana, from 2014 to 2017 were prepared on filter paper. Thirteen years after the removal of drug pressure, chloroquine-sensitive parasite strains with pfcrt K76 have increased nearly to fixation in Begoro, in the forest area (prevalence = 95%), but at a lower rate in Cape Coast, in the coastal region (prevalence = 71%, Z = -3.5, P < 0.001). In addition, pfmdr1 184F-bearing parasites are under strong selection. The pfdhfr/pfdhps quadruple genotype ( IRNG K), associated with SP resistance, is near saturation. Our study identified at a 2 to 10% prevalence pfdhps 581G, which is a sulfadoxine resistance marker that correlates with the failure of SP prophylaxis in pregnancy and which has not been observed in Ghana. The differences in the reexpansion of chloroquine-sensitive strains observed at the two study sites, the stronger SP resistance, and the high prevalence of pfmdr1 184F should be further monitored to inform malaria control strategies in Ghana.

Pooled Deep Sequencing of Drug Resistance Loci from Plasmodium falciparum Parasites across Ethiopia.

Although Ethiopia has an overall lower prevalence of Plasmodium falciparum among countries in Africa, the emergence of drug resistance could seriously hinder elimination efforts. Using samples collected from five therapeutic efficacy studies conducted in 2007-11, we evaluated the prevalence of putative drug resistance mutations in the pfcrt, pfmdr1, and kelch13 genes at the time of those studies, as well as the ama1 gene for genetic relatedness using a pooled amplicon deep sequencing approach. Among all sites, the kelch13 gene showed no mutations, whereas the pfcrt CVIET genotype was fixed in all populations. By contrast, the mdr1 gene demonstrated frequencies of resistant genotypes ranging from 10 to 100% at amino acid position 86 and from 0% to 57.8% at amino acid position 1246. Although we observed a low degree of haplotype sharing between sites, we did observe considerable haplotype sharing within sites over time. This suggests that P. falciparum populations in Ethiopia are isolated and able to persist through time.

Plasmodium falciparum genetic variation of var2csa in the Democratic Republic of the Congo.

BACKGROUND: The Democratic Republic of the Congo (DRC) bears a high burden of malaria, which is exacerbated in pregnant women. The VAR2CSA protein plays a crucial role in pregnancy-associated malaria (PAM), and hence quantifying diversity at the var2csa locus in the DRC is important in understanding the basic epidemiology of PAM, and in developing a robust vaccine against PAM. METHODS: Samples were taken from the 2013-14 Demographic and Health Survey conducted in the DRC, focusing on children under 5 years of age. A short subregion of the var2csa gene was sequenced in 115 spatial clusters, giving country-wide estimates of sequence polymorphism and spatial population structure. RESULTS: Results indicate that var2csa is highly polymorphic, and that diversity is being maintained through balancing selection, however, there is no clear signal of phylogenetic or geographic structure to this diversity. Linear modelling demonstrates that the number of var2csa variants in a cluster correlates directly with cluster prevalence, but not with other epidemiological factors such as urbanicity. CONCLUSIONS: Results suggest that the DRC fits within the global pattern of high var2csa diversity and little genetic differentiation between regions. A broad multivalent VAR2CSA vaccine candidate could benefit from targeting stable regions and common variants to address the substantial genetic diversity.

Matched Placental and Circulating Plasmodium falciparum Parasites are Genetically Homologous at the var2csa ID1-DBL2X Locus by Deep Sequencing.

In pregnancy-associated malaria, infected erythrocytes accumulate in the placenta. It is unclear if in polyclonal infections this results in distinct peripheral and placental parasite populations. We used long amplicon deep sequencing of Plasmodium falciparum var2csa ID1-DBL2X from 15 matched peripheral and placental samples collected at delivery from a high transmission area to determine genetic homology. Despite substantial sequence variation and detecting 23 haplotypes, the matched pairs mostly contained the same genetic variants, with 11 pairs sharing 100% of their variants, whereas others showed some heterogeneity. Thus, at delivery, peripheral and placental parasites appear to intermix and placental genotypes can be inferred through peripheral sampling.

SeekDeep: single-base resolution de novo clustering for amplicon deep sequencing.

PCR amplicon deep sequencing continues to transform the investigation of genetic diversity in viral, bacterial, and eukaryotic populations. In eukaryotic populations such as Plasmodium falciparum infections, it is important to discriminate sequences differing by a single nucleotide polymorphism. In bacterial populations, single-base resolution can provide improved resolution towards species and strains. Here, we introduce the SeekDeep suite built around the qluster algorithm, which is capable of accurately building de novo clusters representing true, biological local haplotypes differing by just a single base. It outperforms current software, particularly at low frequencies and at low input read depths, whether resolving single-base differences or traditional OTUs. SeekDeep is open source and works with all major sequencing technologies, making it broadly useful in a wide variety of applications of amplicon deep sequencing to extract accurate and maximal biologic information.

A deep sequencing approach to estimate Plasmodium falciparum complexity of infection (COI) and explore apical membrane antigen 1 diversity.

BACKGROUND: Humans living in regions with high falciparum malaria transmission intensity harbour multi-strain infections comprised of several genetically distinct malaria haplotypes. The number of distinct malaria parasite haplotypes identified from an infected human host at a given time is referred to as the complexity of infection (COI). In this study, an amplicon-based deep sequencing method targeting the Plasmodium falciparum apical membrane antigen 1 (pfama1) was utilized to (1) investigate the relationship between P. falciparum prevalence and COI, (2) to explore the population genetic structure of P. falciparum parasites from malaria asymptomatic individuals participating in the 2007 Demographic and Health Survey (DHS) in the Democratic Republic of Congo (DRC), and (3) to explore selection pressures on geospatially divergent parasite populations by comparing AMA1 amino acid frequencies in the DRC and Mali. RESULTS: A total of 900 P. falciparum infections across 11 DRC provinces were examined. Deep sequencing of both individuals, for COI analysis, and pools of individuals, to examine population structure, identified 77 unique pfama1 haplotypes. The majority of individual infections (64.5%) contained polyclonal (COI > 1) malaria infections based on the presence of genetically distinct pfama1 haplotypes. A minimal correlation between COI and malaria prevalence as determined by sensitive real-time PCR was identified. Population genetic analyses revealed extensive haplotype diversity, the vast majority of which was shared across the sites. AMA1 amino acid frequencies were similar between parasite populations in the DRC and Mali. CONCLUSIONS: Amplicon-based deep sequencing is a useful tool for the detection of multi-strain infections that can aid in the understanding of antigen heterogeneity of potential malaria vaccine candidates, population genetics of malaria parasites, and factors that influence complex, polyclonal malaria infections. While AMA1 and other diverse markers under balancing selection may perform well for understanding COI, they may offer little geographic or temporal discrimination between parasite populations.

Increased risk of low birth weight in women with placental malaria associated with P. falciparum VAR2CSA clade.

Pregnancy associated malaria (PAM) causes adverse pregnancy and birth outcomes owing to Plasmodium falciparum accumulation in the placenta. Placental accumulation is mediated by P. falciparum protein VAR2CSA, a leading PAM-specific vaccine target. The extent of its antigen diversity and impact on clinical outcomes remain poorly understood. Through amplicon deep-sequencing placental malaria samples from women in Malawi and Benin, we assessed sequence diversity of VAR2CSA's ID1-DBL2x region, containing putative vaccine targets and estimated associations of specific clades with adverse birth outcomes. Overall, var2csa diversity was high and haplotypes subdivided into five clades, the largest two defined by homology to parasites strains, 3D7 or FCR3. Across both cohorts, compared to women infected with only FCR3-like variants, women infected with only 3D7-like variants delivered infants with lower birthweight (difference: -267.99 g; 95% Confidence Interval [CI]: -466.43 g,-69.55 g) and higher odds of low birthweight (<2500 g) (Odds Ratio [OR] 5.41; 95% CI:0.99,29.52) and small-for-gestational-age (OR: 3.65; 95% CI: 1.01,13.38). In two distinct malaria-endemic African settings, parasites harboring 3D7-like variants of VAR2CSA were associated with worse birth outcomes, supporting differential effects of infection with specific parasite strains. The immense diversity coupled with differential clinical effects of this diversity suggest that an effective VAR2CSA-based vaccine may require multivalent activity.

Surveillance for sulfadoxine-pyrimethamine resistant malaria parasites in the Lake and Southern Zones, Tanzania, using pooling and next-generation sequencing.

BACKGROUND: Malaria in pregnancy (MiP) remains a major public health challenge in areas of high malaria transmission. Intermittent preventive treatment in pregnancy (IPTp) with sulfadoxine-pyrimethamine (SP) is recommended to prevent the adverse consequences of MiP. The effectiveness of SP for IPTp may be reduced in areas where the dhps581 mutation (a key marker of high level SP resistance) is found; this mutation was previously reported to be common in the Tanga Region of northern Tanzania, but there are limited data from other areas. The frequency of molecular markers of SP resistance was investigated in malaria parasites from febrile patients at health centres (HC) in seven regions comprising the Lake and Southern Zones of mainland Tanzania as part of the ongoing efforts to generate national-wide data of SP resistance. METHODS: A cross-sectional survey was conducted in the outpatient departments of 14 HCs in seven regions from April to June, 2015. 1750 dried blood spot (DBS) samples were collected (117 to 160 per facility) from consenting patients with positive rapid diagnostic tests for malaria, and no recent (within past 2 months) exposure to SP or related drugs. DNA was extracted from the DBS, pooled by HC, and underwent pooled targeted amplicon deep sequencing to yield estimates of mutated parasite allele frequency at each locus of interest. RESULTS: The dhps540 mutation was common across all 14 sites, ranging from 55 to 98.4% of sequences obtained. Frequency of the dhps581 mutation ranged from 0 to 2.4%, except at Kayanga HC (Kagera Region, Lake Zone) where 24.9% of sequences obtained were mutated. The dhfr164 mutation was detected only at Kanyanga HC (0.06%). CONCLUSION: By pooling DNA extracts, the allele frequency of mutations in 14 sites could be directly determined on a single deep-sequencing run. The dhps540 mutant was very common at all locations. Surprisingly, the dhps581 was common at one health center, but rare in all the others, suggesting that there is geographic micro-heterogeneity in mutant distribution and that accurate surveillance requires inclusion of multiple sites. A better understanding of the effect of the dhps581 mutant on the efficacy of IPTp-SP is needed.

A deep sequencing tool for partitioning clearance rates following antimalarial treatment in polyclonal infections.

BACKGROUND AND OBJECTIVES: Current tools struggle to detect drug-resistant malaria parasites when infections contain multiple parasite clones, which is the norm in high transmission settings in Africa. Our aim was to develop and apply an approach for detecting resistance that overcomes the challenges of polyclonal infections without requiring a genetic marker for resistance. METHODOLOGY: Clinical samples from patients treated with artemisinin combination therapy were collected from Tanzania and Cambodia. By deeply sequencing a hypervariable locus, we quantified the relative abundance of parasite subpopulations (defined by haplotypes of that locus) within infections and revealed evolutionary dynamics during treatment. Slow clearance is a phenotypic, clinical marker of artemisinin resistance; we analyzed variation in clearance rates within infections by fitting parasite clearance curves to subpopulation data. RESULTS: In Tanzania, we found substantial variation in clearance rates within individual patients. Some parasite subpopulations cleared as slowly as resistant parasites observed in Cambodia. We evaluated possible explanations for these data, including resistance to drugs. Assuming slow clearance was a stable phenotype of subpopulations, simulations predicted that modest increases in their frequency could substantially increase time to cure. CONCLUSIONS AND IMPLICATIONS: By characterizing parasite subpopulations within patients, our method can detect rare, slow clearing parasites in vivo whose phenotypic effects would otherwise be masked. Since our approach can be applied to polyclonal infections even when the genetics underlying resistance are unknown, it could aid in monitoring the emergence of artemisinin resistance. Our application to Tanzanian samples uncovers rare subpopulations with worrying phenotypes for closer examination.

Using Amplicon Deep Sequencing to Detect Genetic Signatures of Plasmodium vivax Relapse.

Plasmodium vivax infections often recur due to relapse of hypnozoites from the liver. In malaria-endemic areas, tools to distinguish relapse from reinfection are needed. We applied amplicon deep sequencing to P. vivax isolates from 78 Cambodian volunteers, nearly one-third of whom suffered recurrence at a median of 68 days. Deep sequencing at a highly variable region of the P. vivax merozoite surface protein 1 gene revealed impressive diversity-generating 67 unique haplotypes and detecting on average 3.6 cocirculating parasite clones within individuals, compared to 2.1 clones detected by a combination of 3 microsatellite markers. This diversity enabled a scheme to classify over half of recurrences as probable relapses based on the low probability of reinfection by multiple recurring variants. In areas of high P. vivax diversity, targeted deep sequencing can help detect genetic signatures of relapse, key to evaluating antivivax interventions and achieving a better understanding of relapse-reinfection epidemiology.

Differing patterns of selection and geospatial genetic diversity within two leading Plasmodium vivax candidate vaccine antigens.

Although Plasmodium vivax is a leading cause of malaria around the world, only a handful of vivax antigens are being studied for vaccine development. Here, we investigated genetic signatures of selection and geospatial genetic diversity of two leading vivax vaccine antigens--Plasmodium vivax merozoite surface protein 1 (pvmsp-1) and Plasmodium vivax circumsporozoite protein (pvcsp). Using scalable next-generation sequencing, we deep-sequenced amplicons of the 42 kDa region of pvmsp-1 (n = 44) and the complete gene of pvcsp (n = 47) from Cambodian isolates. These sequences were then compared with global parasite populations obtained from GenBank. Using a combination of statistical and phylogenetic methods to assess for selection and population structure, we found strong evidence of balancing selection in the 42 kDa region of pvmsp-1, which varied significantly over the length of the gene, consistent with immune-mediated selection. In pvcsp, the highly variable central repeat region also showed patterns consistent with immune selection, which were lacking outside the repeat. The patterns of selection seen in both genes differed from their P. falciparum orthologs. In addition, we found that, similar to merozoite antigens from P. falciparum malaria, genetic diversity of pvmsp-1 sequences showed no geographic clustering, while the non-merozoite antigen, pvcsp, showed strong geographic clustering. These findings suggest that while immune selection may act on both vivax vaccine candidate antigens, the geographic distribution of genetic variability differs greatly between these two genes. The selective forces driving this diversification could lead to antigen escape and vaccine failure. Better understanding the geographic distribution of genetic variability in vaccine candidate antigens will be key to designing and implementing efficacious vaccines.